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Past numerical simulations and experiments of turbulence exhibit a liump in the inertial range, 
called the bottleneck effect. In this paper we show that sufficiently large inertial range (four decades) 
is required for an effective energy cascade. We propose that the bottleneck effect is due to the 
insufficient inertial range available in the reported simulations and experiments. To facilitate the 
turbulent energy transfer, the spectrum near Kolmogorov's dissipation wavenumber has a hump. 



I. INTRODUCTION 

i> : 

■ Energy spectrum of turbulent flow is an important quantity. In 1941, Kolmogorov [1] showed that the energy 
\ spectrum E{k) of turbulent flow is 

E{k)^KKon^'^k-'>'^f{k/hi), (1) 

^ ; 

where 11 is the energy flux, Kko is Kolmogorov's constant, kd is Kolmogorov's wavenumber, and the function f{x) 1 
T^lj- I in the inertial range {x ^ 1), and f{x) ^ as x 3> 1. Many experiments and numerical simulations verify this 
(N| . powerlaw apart from a very small intermittency correction. The compensated energy spectrum E{k)k^^^ / Kko is 

flat in the inertial range, and decays in the dissipation range. A careful observation of energy spectrum obtained 
', from recent high-resolution numerical simulations and experiments however show a small hump near Kolmogorov's 
f-^ ' wavenumber kd. The feature is called the bottleneck effect in literature. In this paper we propose an explanation for 

the bottleneck effect. 

I The bottleneck effect has been reported in many numerical simulations and experiments of fluid turbulence. Yeung 
\ and Zhou Gotoh et al. 0], Kaneda et al. Q, and Dobler et al. [H] found the bottleneck effect (hump in 
P5 . the normalized energy spectrum) in their numerical simulations. Saddoughi and Veeravalli @| studied the energy 
spectrum of atmospheric turbulence and reported the bottleneck effect. They observed that the longitudinal spectra 
, have a larger inertial range (around 1.5 decade) but smaller hump, while the transverse spectra have relatively smaller 
^ ■ inertial range (around one decade), but a larger hump. Shen and Warhaft [7], Pak et al. [8], She and Jackson and 
\^ \ other experimental groups also observed the bottleneck effect in fluid turbulence. 

PsJ ' The bottleneck effect has been seen in other forms of turbulence as well. Watanabe and Gotoh and others (lol. [ill, [l^ 
^> \ reported the bottleneck effect in scalar turbulence, Haugen et al. in three-dimensional magnetohydrodynamics, 
. Biskamp et al. [l3| in two-dimensional magnetohydrodynamics, electron-magnetohydrodynamics, and thermal con- 
vection. Lamorgese et al. [lB|, Biskamp et al. [lj|, and Dobler et al. 0| observed that the bottleneck effect became 
more pronounced when hyperviscosity is increased. 

There have been various attempts to explain bottleneck effect. Falkovich [ll| argued that the viscous suppression 
of small-scale modes removes some triads from nonlinear interactions, thus making it less effective, which leads to 
pileup of energy in the inertial interval of scales. Based on turbulent viscosity and the assumption of local energy 
transfer, Falkovich derived the following formula for the correction in Kolmogorov's spectrum: 
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5E{k) = E{k){k/kp)*^^/\n{kp/k), (2) 

. where kp is proportional to the dissipative wavenumber kd, and k <^kp. 

Yakhot and Zakharov [l3| derived energy spectrum using Clebsch variables and showed that the energy spectrum 



IS 



E{k) - KKo^^'^k-'>/^J{k/kd) + Pk-\ (3) 

i.e., correction is of the form k~^ . Theoretical justification for k'^ was argued by Orszag [isj who analyzed the one- 
loop Dyson equation for the propagator G and the velocity correlation function U; the spectrum k'^ was obtained by 
assuming that the response function is dominated by viscous effects. She and Jackson reported an experimental 
result in which they observed the k~^ bottleneck correction; they argued coherent vortex structures to be the reason 
for the bottleneck effect. 



Kurien et al. [l^ extended Kolmogorov's phenomenology to include the effects of helicity. They found a shallower 
k~^^^ energy spectrum at higher wavenumbers by assuming that the helicity transfer time-scales dominate at large 
wavenumbers. In the following discussion we will propose a new mechanism to explain the bottleneck effect. We argue 
that the bottleneck effect is seen when the length of inertial range is insufficient for the energy cascade process. 



II. THE REASON FOR THE BOTTLENECK EFFECT 



The basic idea presented in our paper is as follows: In a fully-developed turbulence, a flux of energy is transferred 
from small wavenumbers to large wavenumbers. This process involves interactions of large number of modes — from 
small wavenumbers to large wavenumbers. The maximum energy transfer from a given wavenumber shell is to its 
nearest neighbour. Still significant amount of energy transfer takes place between somewhat distant wavenumber 
shells (20l. [ill. [2^. [23l. l2i. f25| . Verma et al. [IHl showed using a theoretical arguments that if the inertial-range shells 
are divided in such a way that the mth shell is given by /co(2"/^ : 2('"+i)/4), then in the inertial -range, the normalized 
shell-to-shell energy transfer rates from shell m to shell m+l,m + 2,m + 3 are 18%, 6.7%, and 3.6% respectively. The 
remaining portion of energy flux, which is a huge fraction (^^70%), is transferred to the distant shells. This result is 
in agreement with earlier simulation results (2l| . 

The above arguments imply that for an effective cascade of energy, there must be a large enough range of wavenum- 
bers. Ideally, when Kolmogorov's wavenumber kd 00, Kolmogorov's cascade is setup, and the energy spectrum is 
given by Eq. ([1]). However, if Kolmogorov's wavenumber is not sufficiently large, the cascade process faces difficulty; 
at higher wavenumbers there are not enough number of modes to receive the energy transferred from the smaller 
wavenumbers. To compensate, the wavenumbers near Kolmogorov's scale have a higher energy level. We propose 
this to be the main reason for the bottleneck effect. Note that the energy fed at small wavenumbers fixes the level of 
energy spectrum in the inertial range, and the energy input has to be dissipated at the higher wavenumbers. In the 
following discussion we will present a quantitative arguments to support the above idea. 



Formalism 



The average energy fiux from a wavenumber sphere of radius ko is given by ]2Q, \23, \2L 



n(fco) 



dk / dp{S{k\p\q)) 
k>ko J p<ka 



(4) 



where S(k\p\q) is the ''^ mode-to-mode energy transfer rate" in a triad (p,q,k) with k = p -I- q, and () represents the 
ensemble average. The term S{k\p\q) represents the energy transfer rate from mode p to mode k with mode q acting as 
a mediator. The term {S{k\p\q)) has been computed earlier using standard field-theoretic technique (2^. [27l. f28l. I29I. [30| . 
The procedure to compute n(A:o) is described in the above references, which yields 



n(fco) = k'J^u 



dke r dp r dq^'^ ^^^(^)^('^) + ^^^('^)^(^) + ^^^(^)^(^) 

ko Jo Jlk-pl 



I/*(fc2/3+p2/3 + g2/3) 



(5) 



where v* is the renormalized parameter in the expression of renormalized viscosity |28l . l29l . l30| , and it has been found 
it to be between 0.35 to 0.40. In this paper we take v* — 0.38 [2^. The correlation function C{k) is related to the 
one-dimensional energy spectrum E{k): 



C{k) = 



E{k) 
47rA;2' 



(6) 



and T/s are given in [26 



T3 



kp{xy + 2z^ + 2xyz'^ + z), 
—kp{xy + 2z^ + 2xyz^ + y'^z), 
—kq{xz — 2xy'^z — y'^z), 



where x,y,z are cosines defined as 



p • q = —pqx, q • k = qky, p • k = pkz. 



(7) 



The field-theoretic method mentioned above has certain similarities with the calculations based on the eddy-damped 
quasi-normal Markovian approximation (EDQNM). Both these methods use quasi-normal approximation, and eddy 
or renormalized viscosity. 

In the subsequent subsections we will use the above formalism to compute energy fiuxes using energy spectrum 
obtained from a model and direct numerical simulation. We also estimate the extent of the bottleneck effect using 
energy transfer ideas. 

B. Bottleneck Effect in Energy Spectrum 

We compare our theoretical results with numerical simulation. The simulations have been performed for homoge- 
neous, isotropic turbulence with stochastic forcing at low wavenumbers. These simulations were done at 512'^, 1024'^, 
and 2048"^ grids. Taylor-based Reynolds numbers for these runs were approximately 240, 400 and 700 respectively, (see 
Yeung et al. for details on simulation.) We multiply the numerical energy spectrum with k^^^ /Kko {Kko = 1-58), 
then divide the resultant quantity by its maximum value in the inertial range, and obtain compensated energy spec- 
trum E{k). In the inertial range, E{k) = 1. In Fig. 1 we plot E{k) obtained from direct numerical simulations (DNS) 
done on 512'^, 1024^, and 2048^ grids. A hump appears in all the DNS plots indicating the existence of the bottleneck 
effect in numerical simulations. These results are consistent with earlier numerical results showing the bottleneck 
effect. 

Comparison of the normalized energy spectra for different grid resolutions reveals that the hump is most dominant 
for 512^, and it decreases as the grid size or Reynolds number is increased, a phenomenon observed in earlier numerical 
results as well [1, [13, [HI, IsH . This result indicates that the bottleneck effect decreases with the increase of inertial 
range, thus reinforcing our hypothesis that the bottleneck effect may be due to nonavailability of sufficient range of 
wavenumbers to facilitate energy cascade. Please note that we have quantified the bottleneck effect by the size of the 
hump in the normalized energy spectrum. In individual energy spectrum the size of the hump could depend on the 
energy input rate etc. Also, we observe a hump at low wavenumbers which is due to the forcing at these scales. The 
focus of this paper is on the hump at k ^ kd and we will not analyze the one at the lowest waveumbers. 

Let us compare the above energy spectra with a model energy spectrum for a turbulent fiow ^15 . i3^ . [33| 

E{k) = KKoA{k/kf)U'^/^k~^/^ exp {-ck/kd), (8) 

where 

^s-l-5/3 

^(^) = l + ^.+5/3 (9) 

with forcing wavenumber kj = 2, c = 0.2, and s = 4. Throughout this paper we take Kko = 1-58 (5^, [2§|. Clearly, 
E{k) oc k^ for k < kf, E{k) oc k~^/^ for the intermediate range (fcy < k < kd), and E{k) oc fc~^/'^ exp(— cfc/fc^) for the 
dissipation range (fc > kd). The choice of s = 4 is based on Batchelor's spectra [sJl for smaller wavenumbers. There 
is no hump in the model spectrum because of the choice of its functional form. Here we compare these spectra with 
spectra that show the bottleneck effect in order to see how the latter affects the spectral energy transport. 
Without loss of generality we can take 11 = 1. In Fig. 1 we plot E(k), which is given by 

E{k) = E{k)k'^'^/KKo = A{k/kf ) exp {-ck/kd) (10) 

As expected, E{k) with higher kd produces a larger inertial range. 

In the following subsection we will compute energy flux by substituting the above energy spectra (DNS and model) 
in Eq. |(5|) and compare the results. They provide important clues for the bottleneck effect. 

C. Bottleneck Effect in Energy Flux 

First, we compute the flux Il{k) by substituting the model energy spectrum [Eq. ^] in Eq. ^ with 11 = 1. We 
compute the integral /(fco) (the bracketed term of Eq. [5]) for various values of kd. When kd = oo and A{x) = 1, 
the integral loo = 0.50 independent of ko, implying that the flux is independent of fco for the Kolmogorov energy 
spectrum (E(k) = KKo^^^^k"^^"^). Using loa we flnd the Kolmogorov constant, Kko = 1-58 (this is how Xxowas 
computed in [28\). After this the integral /(fco) is computed using the model spectrum with s — 4, c = 0.2, and 
kd = 100, 1000, 10000. The value of /(fco) starts from at fco = 0, reaches a peak, and then it decays. 



The energy fluxes at various wavenumbers are 

mo) - K'fjiko) (11) 

with Kko = 1-58. Fig. 2 contains plots of n(/co) vs. fco for different values of (fc^, c). The maximum values of n(fco) 
for these cases are Hsted in Table 1. They are all less than 1, but the difference from the actual value (1) is lower for 
larger kd- Theoretically maa;(n(/co)) must be 1 because the energy input at small wavenumber is 1. The reason for 
the decrease in max{Il{ko)) is the lack of modes in the inertial range. This is where the hump in the energy spectrum 
near dissipation wavenumber comes into play. 

After the flux calculation for model spectrum, we compute the flux integral using E{k) obtained from DNS at 512^, 
1024^, and 2048^ grids, and obtain max{IlDNs)- These values are listed in Table 1. The value of max{IlDNs) for 
2048^^ is very close to unity. Clearly the energy spectra obtained from numerical simulations provide a better handle 
on energy flux as compared to the model energy spectrum [Eq. (JH])]. This is because of the higher level of energy 
spectrum (hump) near Kolmogorov's wavenumber in the DNS (see Fig. 1), which makes up for the loss of large 
wavenumber modes. The overall effect is that the energy flux in high- resolution DNS is closer to what is expected in 
an idealized situation when kd oo. Thus consistency with Kolmogorov's theory is achieved. The value max(IlDNs) 
for 512^ is somewhat higher than 1, which may be to approximations made in our theoretical calculations. 

The DNS plots of Fig. 2 are the fluxes computed by substituting the DNS energy spectra in Eq. This exercise 
was done to examine the effects of the bottleneck correction in the flux. In Fig. 3 we plot the normalized energy 
flux computed directly from DNS data on 512^, 1024'^, and 2048^ grids. The two flgures match qualitatively, but not 
quantitatively because of the assumptions made in the fleld-theoretic calculation. The coupHng of wavenumber modes 
in forced, inertial range, and dissipation range is not yet fully understood to be able to resolve n(fc) completely from 
theory (11, [H, [H, [sll • 

In the present subsection we showed that the bottleneck correction near the dissipation range helps in the effective 
transfer of energy flux. In the next subsection we estimate the extent of the bottleneck effect to the mechanism 
proposed in our paper. 

D. Estimation of the Bottleneck Effect 

In this subsection we will attempt to estimate the extent of the bottleneck effect using semiquantitative arguments. 
Because of a lack of complete understanding of the coupling between the forced, inertial, and dissipative scales, this 
is the best we can do at present. 

Verma et al. [IBl and Verma computed the shell-to-shell energy transfer rate from mth wavenumber shell to 
nth wavenumber shell (T™) using 

where Sm,n are the wavenumber range for the mth and nth shell respectively. The wavenumber space is divided into 
various shells logarithmically. In Verma et al. [2^1 and Verma [2^ the mth shell is (2™/* : 2^™"'"^)/'*). 

Verma et al. [25] and Verma [2B] computed T™ in the inertial range using a similar procedure as described in the 
previous subsection. Kolmogorov's spectrum k~^^^ was assumed through out the wavenumber space. They found 
that the energy transfer is maximal to the nearest neighbour, yet signiflcant energy is transferred to other shells. 
For example, the energy transfer rates from mth shell the shells m + 1, m + 2, and m + 3 are 18%, 6.7%, and 3.6% 
respectively. The transfer rate decreases monotonically for more distant shells. 

Let us imagine a wavenumber sphere of radius R somewhat in the middle of the inertial range. Using the shell- 
to-shell energy transfer rates we can compute the energy transfers from the above wavenumber sphere to n shells 
adjacent to the sphere (wavenumber range [R : R* 2"/'']). Simple algebra shows that the above quantity is f2S| 

^^X^mxT^. (13) 

m— 1 

In Table 2 we list Q„ for various values of n. The Table shows that 42.2% of the flux is transferred to the three 
adjacent shells. To transfer 99.1% energy we need 28 shells in the right of the sphere. Therefore, we require large 
number of wavenumber shells for effective energy transfer, and the bottleneck effect is expected if the inertial range 
is insufficient. In this theory, the bottleneck effect would disappear when there are sufficient wavenumber shells to 
enable the complete energy transfer. 



dkk' 



kes„ 



dp 



k—p\ 



-T2C{q)Cik)+T3Cip)C{k) 



'4fc 



(fc2/3 



^2/3 



g2/3) 



(12) 



The energy transfer among the wavenumber shells is antisymmetric, that is T™ = — T^. If we assume the above 
mentioned waveumber sphere to be in the middle of the inertial range, we require approximately 28 x 2 = 56 shells 
for an effective energy transfer. Hence, the inertial range {kmax/kmin) required must be around 2^^/* « 10*. Hence 
our estimate for the minimum length of the inertial range for no bottleneck effect is approximately four decades. The 
range of inertial range in all the experiments and simulations discussed in this paper is less than four decades, and the 
bottleneck effect is observed in all of them. Hence our theoretical estimate is consistent with the present experimental 
and numerical results. We remark that the above estimate of the required inertial range for zero bottleneck effect 
could be an overestimate. A realistic estimate requires a detailed study of energy transfer among modes in the whole 
range: forcing, inertial, and dissipation range. 

After the above estimation of required inertial range to suppress the bottleneck, we move on to estimate the increase 
in the energy spectrum due to the bottleneck effect. Suppose the energy spectrum E{k) till the dissipation wavenumber 
is 

E{k) = i^A'otf /^fc~5/3(^ ^ (^4^ 

where e(fc) = 5E{k) / E(k) is the normalized bottleneck correction. There is a complex interaction between the 
wavenumbers in the forcing, inertial, and dissipation range, which is not yet completely understood. For time being 
we estimate the additional energy transfer due to the bottleneck correction to be of the order oi Tbottieneck ~ H x e{kd). 
Since the energy suppHed at the large-scales has to reach the dissipation scale, and if the number of wavenumber 
shells to the right of above mentioned wavenumber sphere is n, then 

eneck 



Therefore, 

e(fcd) 



n n 



SE{kd) 



1. (15) 

Qn 



1 - 

n 



(16) 



EKol7n{kd) 

Using Zhou (2ll | and Verma et al.'s results that T™ ~ |" — m|^*/'^ for small {n — m), we estimate 

|i « an^^ (17) 

where a is a positive constant. Assuming that the we have equal number of wavenumber shells to the left and right 
of the wavenumber sphere R discussed above, the ratio of Kolmogorov's wavenumber and forcing wavenumbers is 
approximately 

^ - i?e3/4 ^ (2"/4 X T/^), (18) 

fco 

which yields n ~ 1.5\og2 Re, where Re is the Reynolds number based on Kolmogorov's scale. Substituting this 
estimate of n in Eq. (fT6| we obtain 

e(fc<j) - l-a(1.51og2i?e)2/3, (19) 

which is plotted in Fig. 4 for a reference with a — 0.09. The three points represent the [Re,e{kd)) for three DNS 
discussed in the present paper. The choice of a = 0.09 fits best with the DNS values, and it is consistent with our 
Qnequation [Eq. (flTl) ]. The numerical values of DNS fit quite well with the theoretical predictions, however, we 
need more DNS results for a better test of our theoretical estimate of the bottleneck correction. Also, for a = 0.09, 
e[kd) « for n « 37 and Re « 10''. These estimates are in reasonable agreement with our earher estimate of the 
length of the inertial range for zero bottleneck effect. Our prediction of e{kd) is proportional to 1 — const (log i?e)^/'^, 
and it differs from the predictions of earlier theories. 

Please note that the above expression for e{kd) is only a crude approximation, and could be an overestimate. To 
better understand the bottleneck effect we need to understand the coupling among forcing, inertial, and dissipation 
scales, as well as other aspects Hke intermittency. 

The dynamics at the dissipation rate is quite important in the study of the bottleneck effect. This is evident from 
the numerical observations of Lamorgese et al. [l3|, Biskamp et al. [l3|, and Dobler et al. Q who reported that 
hyperviscosity enhances the bottleneck effect. Since the extent of inertial range increases with the introduction of 
hyperviscosity, it may appear that the bottleneck effect should decrease in the presence of hyperviscosity. However, 
that is not the case. This result is possibly because of the shorter dissipation range in the presence of hyperviscosity, 
and the hump in the energy spectrum near kd could help in the inertial-range energy transfer as well as in the 
dissipation of energy. This is an important question to investigate. So far our focus has been on the physics of energy 
transfer in the inertial range. A more detailed study of energy exchange between wavenumbers in the inertial and 
dissipative range is required for a conclusive statement [ssl . l37l |. 



III. CONCLUSIONS 



To summarize, in this paper we investigated the reason for the bottleneck effect in turbulence. The energy is 
suppHed at large scales, and it cascades to smaller scales. Recent numerical and theoretical studies show that even 
though most of the energy from a given wavenumber shell goes to the next wavenumber shell, there is a significant 
energy transfer to the distant wavenumber shells. We showed that an effective transfer of energy flux in the inertial 
range can take place when there is approximately four decades of inertial range. If the inertial range is shorter, a hump 
is created near Kolmogorov's scale (beginning of dissipation range) which compensates for the nonexistence of required 
inertial range. The bottleneck effect is observed in most of the current numerical simulations and experiments. 

The mechanism proposed in the present paper differs from that of Falkovich and Yakhot and Zakharov [l7| . 
Falkovich argued that the bottleneck effect is due to the suppression of nonlinear interactions by dissipative 
modes, and it is present for all dissipative turbulence systems. Falkovich assumes essentially a local energy cascade 
in contrast to both local and nonlocal transfers in our mechanism. In our picture, the energy is transferred to the 
dissipative scales not only from its immediate neighbouring wavenumber shells, but also from the middle of inertial 
range. The energy transfer by Kolmogorov's spectrum requires certain minimum inertial range. If this range is not 
present, the energy levels of the modes near Kolmogorov's scales increase to facilitate the energy transfer. Note 
that if full range of inertial range is present, the last wavenumber shells in the inertial range would transfer only a 
small fraction of energy flux, and there is no bottleneck effect. Our theory suggests that the bottleneck effect will 
disappear if the inertial range is more than approximately four decades. Yakhot and Zakharov and She and 
Jackson [3] obtained bottleneck correction. Our model purely based on energy flux differs from these theories as 
well. Quantitatively, our prediction for the bottleneck correction e{kd) is proportional to 1 — const (log i?e)^/'^, and it 
differs from the predictions of earlier theories. 

Traditional shell models of turbulence assume local energy transfers, and have a large inertial range (2^5-20 ^ 
lQ5-ioy rpj^g bottleneck effect is generally not observed in the shell models. However Biferale and kerr [38] report 
the bottleneck effect in a shell-model (n = 15) based on Kerr-siggia model. So shell models with small inertial range 
could show bottleneck effect, but the bottleneck effect in shell model is in the spirit of Falkovich's mechanism; there 
is not enough dissipative scale to dissipate the cascaded energy. 

The "real" turbulence however involves local as well as nonlocal energy transfers that are not simulated in local 
shell models of turbulence. The recent nonlocal shell models [s^ attempt to model these features of turbulence, and 
it will be interesting to investigate bottleneck effect in the nonlocal shell models. We remark that the field-theoretic 
calculation presented in this paper is more fundamental than the shell model, and some of it's features are same as 
the shell model. Still it is instructive to independently investigate bottleneck effect using nonlocal shell model. 

Many important and unresolved issues are involved in the study of the bottleneck effect. We need to fully understand 
the nonhnear coupling between the forcing inertial, and dissipative range (see Alexakis et al. [s^, DebHquy et al. 
[3^ . Verma et al. [1^, Brasseur and Wei (321 for some of the recent attempts). The vortex interactions, intermittency 
etc. also come up in the study of the bottleneck effect, and we need to understand them better as well. 

The energy transfer in the turbulence of passive scalar and magnetohydrodynamics follows similar patterns as in 
fiuid turbulence. The energy transfer is forward and local, yet significant range of inertial-range is required for effective 
energy transfer ji^l. Hence we expect the bottleneck effect to be present in these systems as well. These projections are 
consistent with a strong bottleneck effect observed in numerical simulations of Watanabe and Gotoh [10| and Yeung 
et al. (ill . [1^ for passive-scalar turbulence, and those of Haugen et al. [l3| for MHD turbulence. The bottleneck 
effect has been observed in electron-magnetohydrodynamic (EMHD) turbulence, and two-dimensional turbulence (see 
Biskamp et al. and references therein), but its cause is possibly more complex. It has been observed that the 
bottleneck effect along the transverse and longitudinal directions are different |6(|; this result still lacks satisfactory 
explanation. Future developments in theoretical turbulence will possibly resolve some of these issues. 
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Table I: The maximum values of energy fluxes for the model energy spectra [Eq. ^] with c = 0.2 and kd = 100, 1000, 10000, 
and for energy spectra obtained from numerical simulations on 512^,1024^ and 2048^ grids. For Kolmogorov's spectrum 
max(n{ko)) = 1. 
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Table II: The energy transfer rates from a wavenumber sphere in the inertial range to n shells adjacent to the sphere (Qn/II = 
J2m=i * fo"" various n's. 
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Figure 1: The normalized compensated energy spectra E{k) — E{k)k^^^ / Kko vs. k for a model energy spectrum [Eq. 
with c = 0.2 and kd = 100, 1000, 10000, and from numerical simulations on 512^, 1024^ and 2048^ grids at steady state. We 
take the energy flux 11 = 1 in the inertial range, so that E{k) = 1 in the inertial range. 




Figure 2: The computed flux n(A;o) using Eq. © for a model energy spectrum [Eq. ©l with c = 0.2 and kd = 100, 1000, 10000, 
and from numerical simulations on 512"^, 1024"^ and 2048^ grids at steady state. 
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Figure 4: A plot of our estimated normalized bottleneck correction e(fcd) — 5E{kd)/E{kd) [Eq. lfT9|l ] as a function of Reynold 
number Re. The three points (*) represent [Re,e{kd)\ for DNS on 512^, 1024^ and 2048^ grids at steady state. 



